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We propose a scheme that allows to analytically determine the three elementary cross sections 
and connect the solutions to the relative sign between the proton and the neutron spin scattering 
amplitudes once the measurements of total event rate from three appropriate targets become avail- 
able. In this way it is thus possible to extract the maximum information on the supersymmetric 
parameter space obtainable from direct detection experiments, in the case that the dark matter 
particle is the lightest neutralino. Our scheme is based on suitably normalized form of the isospin 
momentum dependent structure functions entering in the spin-dependent elastic neutralino-nucleus 
cross section. We compare these functions with the commonly used ones and discuss their advan- 
tages: in particular, these allow in the spin-dependent cross section to factorize the particle physics 
degrees of freedom from the momentum transfer dependent nuclear structure functions as it happens 
in the spin-independent cross section with the nuclear form factor. 

PACS numbers: 95.35.+d, 12.60.Jv 



43 

Or 
i ' 

<L> : 
43 . 



> 
oo 
O 



X 



I. INTRODUCTION 

Dark matter (DM) is the dominant matter component 
in the Universe. Its nature, however, remains elusive. 
To settle this point is an issue central to both particle 
physics and astrophysics. If the dark matter halo sur- 
rounding our galaxy is composed of weakly interacting 
massive particles (WIMP), one expects them to interact 
at some level with terrestrial detectors depositing a 
tiny amount of energy on the recoiling nucleus. From 
such experiments one hopes to extract the elementary 
WIMP-nucleon cross sections which should provide infor- 
mations on the nature of the DM particle. Many exper- 
iments, currently running or planned, aim at achieving 
this goal. 

One of the most studied WIMP candidate is the light- 
est neutralino in the minimal supersymmetric standard 
model (MSSM) with imparity conservation [2tJ|- Neu- 
tralino scattering with nucleons [BHl] arises at the ele- 
mentary level from its interactions with quarks. Dia- 
grams with the exchange of scalar particles, neutral CP- 
even Higgs bosons and squarks, give a scalar interac- 
tion that in the non-relativistic limit results into a spin- 
independent (SI) interaction. The WIMP-proton and the 
WIMP-neutron SI cross sections are almost equal, thus 
the scattering with the nucleus is coherent, i.e. the SI 
cross section is proportional to A 2 , A being the atomic 
mass number. Furthermore, the interaction mediated by 
the Z boson and squarks induces a spin-dependent (SD) 
WIMP-nucleon scattering and a SD interaction with the 
whole nucleus. 

Assuming that the above nucleon cross sections are of 
the same order, the SD event rate is expected to be much 
smaller than the SI rate, even in the case of fairly light nu- 
clei, due to the mentioned A 2 dependence of the SI cross 
section @ . Given that the SD cross sections carry funda- 
mental information regarding the the underlying theory, 



the natural question is, then, whether direct detection 
experiments can measure them. 

The fact that the above discussed elementary quan- 
tities are all that can be constrained or extracted from 
direct detection experiments using three types of nuclear 
targets for any fixed WIMP mass was argued in Ref.fToj 
(see also Ref. [n|). 

In this paper we propose a scheme that allows to an- 
alytically determine the three elementary cross sections 
and connect the solutions to the relative sign between the 
proton and the neutron spin scattering amplitudes once 
the measurements of total event rate from three appropri- 
ate targets become available, thus extracting the maxi- 
mum information on the SUSY parameter space available 
from direct detection experiments. 

Other multi-couplings/multi-targets approaches and 
analysis are found in literature [ll| . In Ref. [l2[ a scheme 
for an analysis of the DAMA results considering both the 
SI-SD cross sections was proposed. In Ref. [13| the au- 
thors developed a procedure for deducing upper limits 
for both the spin dependent couplings. An analysis of 
the experiments with SD sensitivity, in terms of both 
the SD proton and neutron cross section, was performed 
for example in Refs. [34I3- in some recent pap ers 
the Bayesian and frequentist statistical approach [I8l.ll9j 
have been applied to the multi-target analysis to esti- 
mate the uncertainties in the measurements of the basic 
quantities. 

Our scheme is based on a suitable normalized form of 
the momentum dependent spin structure functions (SSF) 
entering in the SD nutralino-nucleus cross section which 
differs from the standard one. Though this functions 
were used by one of the present authors in previous pub- 
lications |20l - [27| . in this work, we explicitly compare the 
two formalisms and emphasize the advantages of the nor- 
malized SSF. 

The plan of the paper is as follows. In Section|TT]we dis- 
cuss our scheme and the main results. In Section lllll we 
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compare with the predictions of the constrained MSSM 
(CMSSM). Further discussions and comments are given 
in Section Qvl Finally, in the Appendices we discuss and 
demonstrate the advantages of the normalized spin struc- 
ture functions over the standard formalism. 



II. AN ANALYTICAL SCHEME FOR 
EXTRACTING THE ELEMENTARY CROSS 
SECTIONS 

For a nucleus with mass number A the total event rate 
per unit of mass per unit of time is given by [27] 



given by o p ® = 3(^ p 2 /7r)|a P! „| 2 . Using this expression^ 
and introducing g = ±1, we get 
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Here pi — 0.3 GeV/cm 3 is the local dark matter den- 
sity density in the halo, m x the WIMP mass, m p w m n 

the nucleon mass, i.e. Am p is the nuclear mass, of^ (0) not va lid using the commonly used functions 5y j| 

and ofR (0) are the total SI and SD WIMP-nucleus cross 



The factors t SI ' SD are the convolution of the nuclear 
form factor and of the momentum transfer dependent 
SSF with the WIMP velocity distribution, assumed to be 
the standard truncated Maxwellian ve locit y dist ribu tion. 
The root mean squared velocity is \J (v 2 ) = / 2vq = 
280 km/s with v = 229 km/s the velocity of the Sun 
around the center of the galaxy. We neglect the small 
time dependent effects (modulation) associated with the 
motion of the Earth. 

It is worth to remark here that the factor t$D in the 
SD part of the total rate in Eq. (JTJ) does not depend 
on the particle physics couplings a PiTl but only on the 
momentum transfer dependent SSF. This factorization is 
achieved using the SSF defined in Refs. [20l - [27l | while 



sections in the zero momentum transfer limit (ZMTL). 

Assuming the near equality of the proton and neutron 
SI cross sections as it is the case for neutralino, in terms 
of the elementary WIMP-nucleon cross section a SI , the 
SI WIMP-nucleus cross section reads 
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with Hp (ha) the WIMP-proton (nucleus) reduced mass. 

Moving to the SD case, the ZMTL spin nuclear matrix 
elements are defined as 



|48J . We discuss this issue in full details in Appendix [S] 
In view of the fact that the spin structure functions 
have been obtained in the context of the nuclear shell 
model, we deviate from the usual practice of dark matter 
calculations [28| and use also for the SI cross section a 
form factor obtained in the context of the shell-model. 
Details on t SI ' SD are given in Appendix IB! 

We then proceed by first rewriting Eq. (fTJ) in a form 
where only quantities with the dimension of a cross sec- 
tion appear. To this end we define the factor 
1 
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with J the total angular momentum of the nucleus in 
the ground state and (S p>n ) the expectation values of the 
spin of the proton and neutron groups. The ZMTL SD 
WIMP-nucleus cross section can be written in the proton- 
neutron representation using Eqs. ([3]), (IA3|) . (|A4[) . (|A5I) . 
The SD WIMP-nucleon scattering amplitudes a p and a n 
in general can have opposite sign (see Section ITVT) . thus 
we write Q 



where <7q = 10~ 8 pb is the current experimental sen- 
sitivity to the SI cross section. Rq is thus fixed to 
Rq ~ 1.6 x 10" 5 kg" 1 yr _1 . Furthermore we define the 
abbreviations 
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For the scope of this paper it is convenient to rewrite 
Eq. Q in terms of the elementary cross sections that are 



1 Depending on the nucleus, also the spin matrix elements can 
have opposite sign. The relative sign thus would appear between 
|f2p(2p| and [f2^a n |. Anyway, the sign of £lp n is known from 
nuclear physics calculations. In particular, in case of 127 I, 73 Ge, 
19 F used below, the relevant matrix elements are all positive. 



2 Different normalizations for the SD cross sections can be found in 
literature with the scattering amplitudes still indicated as Clp n ■ 
For example, in Ref. Q| one finds o"p^ = 12(^ p 2 /7r)|a Pin | 2 while 
in Ref. [Tj| = 24G|,(/jp 2 /7r) |a Pi „ | 2 ; in these cases other 

factors appear in Eq. (|4j . For our scope it is useful to have com- 
pact expressions, thus our |o Pi7 i| 2 include all numerical factors 
except the factor 3, that in last analysis comes from the spin 
average. This normalisation seems to be same adopted in the 
code DarkSusy used in Section llVI for numerical calculations. 
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Leaving the dependence on the WIMP mass understood, 
Eq. ([1} can be rewritten in the desired form: 



S A = o 



si 



A l a SD 



(11) 



We now consider three nuclei with the following char- 
acteristics: A\ for which both proton and neutron spin 
matrix elements are relevant, a proton spin favouring tar- 
get Ai with dp 2 3> f2^ 2 and a neutron spin favouring 
target A3 with f2^ 3 3> . A system of nuclei matching 
these conditions is given by three odd mass nuclei cur- 
rently employed in many experiments: 127 I as the target 
Ai, 19 F as the A 2 nucleus and 73 Ge as the A3. 

Detailed nuclear models calculations for the spin ma- 
trix elements (Sp,n) and form factors have been carried 
out in the past years. For 19 F we use the spin matrix 
elements of Ref. glj] that give fl™ = 1.646 and fl™ = 
-0.030; for 73 Ge the results of Ref. ft 73 = 0.066 and 
fi 73 = 0.836; finally for 127 I the set obtained with the po- 
tential Bonn-A from Ref. [13 that furnishes ti 127 = 0.731 
and il} 27 — 0.177. The hypothesis on the spin matrix el- 
ements are thus well satisfied. 

Neglecting the neutron contribution in Ai and the pro- 
ton contribution in A3, from Eq. ([lip we have the system 
of equations: 



Though the solutions do not depend on g explicitly, 
for each sign the solution is unique. In fact, each of the 
three equations of the system in Eq. (TT2"j) is the sum of 
a SI with a term that depends on the SD cross sections, 
say a SI + fi, i = 1,2,3. Obviously, for each solutions 
set, the three terms fi have to be equal, that is, it must 
be fi — fi — f$. However, the square in f\ in the first 
equation contains an interference term that can be posi- 
tive or negative because depends on g. This means that, 
for example, the set s+ cannot satisfy at the same time 
fi(e = +1) = fi = fa and h(g = -1) = fi = f 3 . 
Given that the same argument applies to s_, each set in 
Eq. ([20]) is thus associated to one sign of g. 

In the specific case of 127 I, 73 Ge, 19 F, the factors t SI ' SD 
were calculated for E m in — and neutralino masses up 
to 200 GeV in Ref. [27[. Using those results we find that 
the parameters defined in Eq. (fl~3"|) are very similar for 
all the masses, a ~ 4 x 10~ 3 and b ~ 1 x 10~ 2 , that is 
s« 1 and tel. Considering the cases 1S1 ^ Si ~ S3, 
Si 3> (Si, S3) and Si <C (Si, S3), one easily finds from 
Eqs. flTll and QU) that A ~ 0, A ~ 2jS 2 , A - 27<S 2 <S 3 , 
respectively, thus the discriminant is always positive. 

Some of the cross sections in s + and s_ may turn out 
to be negative for some values of the rates, thus non- 
physical. We have varied numerically the three rates in 
four orders of magnitude and in Fig. Q] we show the be- 
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Solving for a~ and from the second and third 

equations and substituting them into the squared of the 
first, we obtains a second order equation in a SI , namely 
^4(cr S7 ) 2 - 2B(a SI ) +C = 0. Introducing the parameters 
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and the abbreviations 

a = 1 — a — b, 
/? = S Al - aS A2 
7 = 2ab, 

the coefficients are found to be 

A = a 2 - 2 7 , 
B = a/3 - 7 (S A2 
C = /3 2 - 2 7 5 A2 5 A3 . 
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The two sets of solutions, indicated as s+ and s_ respec- 
tively, are thus: 
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Figure 1. Correlated plane (R( 127 l), R( 19 F)) (left column) 
and (i?( 127 I), i?( 73 Ge)) (right column), showing the behaviour 
of the signs in Eq. (SB) for m x = 30 GeV (top) and 200 
GeV (bottom). For example, fixing a value of _R( 127 I) in the 
abscissas of the two columns, one can read at the same time 
on the ordinates, the allowed ranges of values for i?( 19 F) and 
i?( 73 Ge). The cross sections of s+ are all positive for values 
of the rates in the coloured areas while the values of s_ are 
all positive only in the restricted central turquoise and orange 
areas. In the white non- physical areas both s+ and s_ have 
at least one negative value. 
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Figure 2. The solutions, Eq. ([20]|, of the system Eq. (fT2|) for 
the nuclei 127 I, 73 Ge, 19 F as a function of i?( 19 F): s+, solid 
lines, and s_, dashed lines. We fix R( 127 l) = R( 73 Ge) = 1 
kg -1 yr -1 and three values of the neutralino mass. 



haviour of the signs of Eqs. (|20|) in the correlated planes 
(i?( 127 I),i?( 19 F)), left column, and (i?( 127 I), i?( 73 Ge)), 
right column. In the central turquoise and orange re- 
gions around the diagonal, all the three cross sections of 
,s + and all the three cross sections of s_ are positive. In 
the grey and yellow areas, s+ have positive values but s_ 
can have at least one negative cross section. In the white 
areas both the sets have at least one negative cross sec- 
tion. As an example, suppose that a rate i?( 127 I) = lfT 1 
kg -1 yr~ x is measured for the case m x = 30 GeV. In 
the interval of values 10 -3 < i?( 19 F) < 6 kg -1 yr~ x and 
10~ 3 < i?( 73 Ge) < 10 kg -1 yr -1 the solutions s + are 
positive but the solutions s_ are all positive only in the 
restricted ranges 10" 3 < i?( 19 F) < 8 x 10" 1 kg" 1 yr -1 
and 10" 2 < i?( 73 Ge) < 4 kg" 1 yr" 1 . The variations 
of the allowed regions passing from m x = 30 GeV to 
m x = 200 GeV are small. 

We plot in Fig. [2] the extracted cross sections, Eq. (|20l) . 
as a function of _R( 19 F) fixing the rates in 73 Ge and 127 I 
near the experimental limits, _R( 73 Ge) = i?( 127 I) = 1 
kg -1 yr -1 and three values of the neutralino mass. Note 
that the abscissa range of the curves corresponds to the 
interval of values of i?( 19 F), where the both sets of solu- 
tions are positive, that can be read from the left column 
of Fig. [2] for i?( 127 I) = 1 kg" 1 yr" 1 . The two solutions 
for cTp D are almost overlapping: this is because we plot 
the cross sections as a function of the rate in 19 F that is 
the nucleus more sensitive to this cross section. For the 
same reason a^ D changes by orders of magnitude with 



i?( 19 F) while the variations of a SI and a^ D are less dra- 
matic. We find numerically that, for the chosen nuclei in 
the considered range of rates in Fig. O s + ( s_) is the 
set of solutions associated to g — —1 (g — +1). 



III. COMPARISON WITH THE CMSSM 

In the CMSSM [U all the soft breaking terms at the 
weak scale are obtained by evolving the renormalization 
group equations of the model from the common scalar 
mass mo, the common gaugino mass mi/2, the common 
trilinear scalar coupling Aq, assigned at the gauge unifi- 
cation scale. The running further depends on the ratio 
of the vacuum expectation values of the Higgs doublet, 
tan /3, and the sign of Higgs mixing term [i, which we 
take to be positive. 

Assuming that the neutralino constitutes all the DM 
relic density as deduced from WMAP 34 1 measurements, 



0.094 < nh 2 < 0.128 at 3a, the CMSSM parameter 
space is heavily restricted. We further impose (32l . [33[ 
the LEP bound on the light Higgs mass m/, > 114 GeV, 
and other phenomenological constraints using the code 
DarkSusy [35|]. The scan of the parameter space is done 
fixing Aq = 0, tan/3 = 10, 50 and varying the scalar 
mass and gaugino mass up to the values mo < 4000 GeV, 
mi/a < 1500 GeV. 

In Fig. [3] we plot the SI neutralino-proton cross section 
(the neutron SI cross section may differ from that for 
the proton by less than a few percent), the proton and 
neutron SD cross sections and their ratio as a function 
of the neutralino mass. In all the panels the points on 
the bottom arise from points of the parameter space in 
the stau-coannihilation region where the lightest stau has 
mass nearly degenerate with the lightest neutralino. The 
points on the top are in the focus point regions of the pa- 
rameter space, large mo, where the neutralino has a rele- 
vant higgsino component, and some points, at large tan /3, 
in the funnel region, where the neutralino mass is al- 
most half of the mass of the neutral pseudo-scalar Higgs. 
We also plot in the first panel the upper limit curve on 
the SI cross section obtained by the XENON100 experi- 
ment [3^ | . which, at present, has provided the strongest 
constrain on a si- The points with the highest cross sec- 
tions at low neutralino mass, mostly in the focus point 
region, are already excluded. 

The proton and neutron SD scattering amplitudes, 
a p = Ylqdq^l an d a n — ^2 q 9qAq^ l \ are determined 
by the low energy effective couplings g q and Aq , the 
spin fraction of the nucleon carried by the light quarks 
q = u,d, s. We adopted for the latter quantities the 
central values obtained by the HERMES collaboration 
for the proton [37], Au& = 0.842, Ad& = -0.427, 
As^ = -0.085 and zW") = AS*>) , Ad (n ^ = Au^\ 
/\ s ( n ) = As^ from isospin symmetry for the neutron. 
The ratio /cr„ is almost constant and it is greater 
than one in the focus point and funnel strip and less than 
one in the stau-coannihilation regions. 
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Figure 3. Neutralino-nucleon cross sections as a function of 
the neutralino mass in the constrained minimal supersym- 
metric standard model with Aq — 0, n > 0, mo < 4000 GeV, 
m 1/2 < 1500 GeV, tan/? = 10 and tan/3 = 50. The points 
satisfy the WMAP 3a bound 0.094 < Qh 2 < 0.128 on the 
relic density. The blue points are for tan/3 = 10, the turquoise 
ones for tan/3 = 50. The red curve is the XENON100 upper 
limit on the SI cross section. 



proton and neutron spin dependent, and that the solution 
are connected with the relative sign between the proton 
and neutron spin dependent scattering amplitudes. 

These three nuclei were chosen for several reasons: the 
nuclear spin matrix elements are known with sufficient 
degree of precision and their values allow to write the 
equations for the total rates in the form of Eq. ([T2"|) ; the 
parameters defined in Eq. (fT5|) are both <C 1 and this 
ensures that the solutions given by Eq. (f2Tj)) are real for 
all values of the measured rates; finally, they are going 
to be employed in experiments with much more massive 
targets. For example, the COUPP experiment, bubble 
chamber based on CFrI I39j |. which has a proposal for 
a ton scale experiment [4fJ] cannot measure the energy 
spectra but only counting rates (see also the discussion 
in Ref. [IH). This experiment contains the nuclei I and 
F employed in our analysis. If, in the future, precise 
rates are measured by this experiment, in conjunction 
with a ton scale experiment employing Ge like Super- 
CDMS/GEODM 03 and EURECA g| projects, our 
results will directly allow the extraction of all the cross 
sections, even if the recoil energy spectra are not avail- 
able. 

Certainly, due to the great variety of nuclei employed 
in present and planned experiments, other combinations 
of target nuclei can be used, depending on the availability 
of reliable data. If the nucleus A\ has sensitivity only to 
the SI interaction, it is a = b = 0. From Eqs. flt p -(fig] l it 
is found A=l,B = S Al ,C = B 2 thus A = B 2 - AC = 
and the two sets in Eq. (|20p reduce to the unique solution: 



Moreover we find that a p is always negative and a n al- 
ways positive and this crucially depends on the values of 
the spin structure functions of the proton [38]. The fun- 
damental cross sections are thus provided by the set s + . 
Comparing the predictions of s+ in Fig. [2] solid lines, for 
m x = 200 GeV (in the CMSSM masses below - 50 GeV 
are excluded) with Fig. [3] we see that the extracted cross 
sections are compatible with the cross sections predicted 
by the CMSSM in the stau coannihilation region for a 
total rate in fluorine at the level of 5 x 10~ 2 kg -1 yr _1 . 



IV. DISCUSSION 

Other realizations of the supersymmetry breaking or 
supersymmetric models other than the CMSSM may pre- 
dict different relations between the elementary cross sec- 
tions; it is however clear that the SD interactions carry 
a lot of information on the fundamental theory. In spite 
of this fact, the SD event rate can be small also in light 
nuclei thus rendering the measurement of the SD WIMP- 
nucleon cross sections difficult. 

We have shown that, if the total event rates are mea- 
sured in 127 I, 73 Ge and 19 F targets, the system given by 
Eq. (fTTj) allows the extraction of all elementary nucleon 
cross sections, i.e. the spin independent as well as the 



-SI _ <; ^SD _ Sa 2 - OAi SD _ ^As ~ gAi 

(21) 

These are the expected formulas for the simpler case 
where A\ is an even-even nucleus, like Ar employed in 
large mass experiments like WArP [4j|, ArDM (JSJ and 
the future project DARWIN [45|. In this case it is pos- 
sible to extract the SI cross section cleanly but the con- 
nection of the solutions with the sign is lost. However, it 
can be recovered using 127 I with other combinations like 
I, Ar, Ge or I, Ar, F or Xe instead of Ge. 

One last comment regarding the normalized SSF dis- 
cussed in Appendix[X]is in order. In the usual formalism 
of the SD rate employing Eqs. (|A1[) . (|A2j) . in general it 
is not possible to convert an experimental upper limit 
to an upper_pound on ct^d independently of neutralino 
properties (4(| because the nuclear momentum dependent 
degrees of freedom are not decoupled from the particle 
physics couplings. As a matter of fact, has become a stan- 
dard procedure to use the method of Ref. [l3| to set limits 
on the cross sections in a WIMP-independent way. This 
independence is anyway achieved using the ZMTL total 
cross section in Eq. (|A5|) and not the full g-dependent 
cross section: this can be a bad approximation in the 
case of heavy nuclei and/or heavy neutralinos. The cor- 
rect q-dependence, at least for the nuclei for which full 
nuclear physics calculations are available, can be incor- 
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porated in this method employing the normalized SSF 
which allow the factorization of the particle physics de- 
grees of freedom from the g-dependent nuclear physics 
structure functions. 
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Appendix A: The normalized spin structure 
functions 

The calculation of the elastic neutralino-nucleus full 
momentum transfer dependent cross section is conve- 
niently performed in the isospin basis. In the formalism 
introduced by in Ref. [47] and reviewed in |48|, which 
is the commonly used in the particle physics and astro- 
physics literature (see the reviews 0t3|)i the differential 
cross section is written as 



dq 2 



1 



4(/M«) ; 



S(q) 

s(oy 



(Al) 



The momentum transfer dependent SSF enter in the fac- 
tor 



5(g) = alSoo(q) + aoaiS i(q) + a?5n(g). 

In the ZMTL reduces to 
2J+ 1 



(A2) 
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7T 

ao{{S p ) 



J(J + 1) 

+ (S n )) + ai((S P ) - (S n )) 
23 



-i 2 
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being a^i the isoscalar and isovector coupling determined 
by the particle physics model and related to the proton- 
neutron representation by the relations 



«o — a p + On, ai — a p — a n . 
In the ZMTL the total cross section is 



2J + 1 



5(0). 



(A4) 



(A5) 



The single functions 5^(g), i,j = 0,1 denoting the 
isospin channels, that are obtained from nuclear shell- 
model calculations, are not normalised to one, but the 
normalisation of the "spin form factor" is achieved by 
the ratio 5(g)/5(0). 

In the formalism of Refs. [13, l2lj |. the SSF are given 
by: 



Fij(u) 



E 

A,fe 



(A6) 



where f^'^ (u) are the matrix elements of the multi- 
pole expansion of the nuclear spin operators evaluated in 
a shell-model multiparticle basis and S7^°' 1 ' ) (0) = Oj are 
the ZMTL values. Details can be found in Ref. [2lj . 

The dimensionless variable u is given by u = q 2 b 2 /2 
where b is the nuclear oscillator size parameter that de- 
pends on the atomic mass number of the particular nu- 
cleus, b = 1 fm A 1 / 6 w 5 x 10~ 3 MeV -1 A 1 ' 6 . The differ- 
ential neutralino-nucleus cross section in the laboratory 
frame reads 



d<r SD 

(A) 



1 



du 2(fi A bv) 2 2J + 1 



4/4 ty \ 



(A7) 



with 



2 J + 1 

-{al^l 2 Q F m {u) + 20001^0^1^01 (u) 



16tt 



+a?n?Fn(u)). 



(A8) 



The SSF defined Eq. (|A6|) are normalized to one at u = 
because the ZMTL values f2o,i has been factored out thus 
explicitly appearing in Eq. (|A8|) . These are related to the 
fl p , n defined in Eq. ([3]) by 



fig — -(- ^1 — ^p — 



(A9) 



The formulas connecting the two sets of SSF are: 
2 7-1-1 

s oo(q) = -r^-ngJboM, (Aio) 

(All) 
(A12) 



2J + 1 

5oi(g) = — ^o^iFoi(u), 

2-7 +W, 



Su(q) = -r -j— n^ u (u). 
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We are in the position to appreciate the differences 
between the two formalisms. In the case of the light 
nucleus 19 F, the functions Ffj are given in an analytical 
form in Refs. [H M as 



f£ 9f) (u) = e~ v > f' ;; ,» 



My^ f (fc)„,fc 

k=Q 



(A13) 



For the isotope 73 Ge, the authors of Ref. [2!| furnish 
a polynomial fit to the calculated SSF of the following 
form: 



fc=0 



(A14) 



while for 127 I, the authors of Ref. [30( fit their functions 
(here we use the ones obtained with the Bonn-A nucleon- 
nucleon potential) with 



k=0 



(A15) 



in terms of the variable y = (go/2) 2 = u/2. The coeffi- 
cients of the polynomials are found in the cited papers 
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Figure 4. The momentum transfer dependent spin structure 
functions Sij (left column) and Fij (right column), from top 
to bottom for the nuclei 19 F, 73 Ge, 127 I. 



Given that F o(u) ~ F i(u) ~ Fu(u) we can describe 
the momentum dependence of the cross section in terms 
of only one structure function, for example the one as- 
sociated with the isovector channel Fn(u). Considering 
the intrinsic uncertainties in the nuclear physics calcu- 
lations, in the WIMP velocity distribution function and 
in the particle physics couplings, this is by far a well 
motivated approximation with negligible error in the cal- 
culation of the cross section. We thus rewrite Eq. (|A8[) 
as 

2 J + 1 

T(u) = ——{a n + a^fFniu), (A17) 

l07T 

and going to the proton-neutron representation using 
Eqs. dM}, CMl, (IMJ), CM), we finally find 



dcr SD 



1 



(A18) 



Eq. (IA18[) is equivalent to, but simpler than, Eq. (|A1I) . 
We summarize and highlight the advantages of the nor- 
malized SSF: 

1. the three F^ are normalized to one at zero mo- 
mentum transfer and are always positive definite, 
while the interference term Soi, depending on the 
nucleus, can be negative (see the case of 73 Ge); 

2. the three F^ are to a very good approximation 
identical, thus one has to deal only with one SSF 
instead of three; 

3. in Eq. (|A"7|) . differently from Eq. (|AT|) . the particle 
physics and nuclear physics momentum dependent 
degrees of freedom are thus completely factorized. 



and in the review 0]. Using Eqs. (|Al0ll - (|Al5|) and the 
ZMTL Sl ,i obtained with Eqs. © and ([A"9]l. we plot in 
Figure @] the two sets for the nuclei 19 F, 73 Ge, 127 I as 
a function of the variable y. Introducing the quantity 
Qo = (^Amp)^ 1 the variables u and y are related to the 
recoil energy through the relation 

Er = uQ = 2yQ = 2y x 40 /T 4/3 MeV. (A16) 

As it is evident, while the Sij look quite different, the 
F^ are practically identical in the interesting range < 
y < 1: using Eq. (IA16|) . y = 1 corresponds to Er = 1578 
keV in 19 F, E R = 262 keV in 73 Ge and E R = 125 keV in 
127 I. We remark that while for the light nucleus 19 F the 
functions in Eq. (|A13[) are an exact result of the shell- 
model calculations, the functions in Eqs. (|A14[) . (|AI5[) 
for the more complex nuclei 73 Ge, 127 I, are fits to the 
results of the shell-model calculations which cannot be 
cast in a simple analytical form. Furthermore, in a more 
recent large-scale shell-model computation performed in 
Ref. [4!| the authors calculate the Fij{u) SSF for the 
heavy nuclei 127 I, 129 Xe, 131 Xe, 133 Cs and find the same 
behaviour (see Fig. 3 and Fig. 4 of [49]). 



Appendix B: The factors t bI,bn 

Starting from Eq. (fA7]l the t SD factor in Eq. (TJJ is 
defined as 

,-=7'™ (B1) 

«mi» !) mj „{u) 

As stated in Section HH we neglect in this work the 
small effects associated to motion of the Earth and con- 
sider the standard WIMP Maxwellian velocity distribu- 
tion truncated at the escape velocity v esc *5Q| . The inte- 
gration limits are: v ma x = v esc , v min (u) = v /(u/(2^& 2 ), 
u m ax = 2(^ A bv esc ) 2 and u min = E min /Q Q , where E min 
is set by the energy threshold of the detector. 

The factor t SI is given by the same formula replacing 
Fu(u) with |F(m)| 2 , where F(u) is the nuclear form fac- 
tor [20, HU, HH . For consistency we use the same nuclear 
model, i.e. the same orbitals, the same interaction and 
the same size parameter, in obtaining the spin structure 
functions and the coherent form factor. In the context 
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of the shell-model, the coherent form factor arises: (1) 
from the contribution of the A caTC = A — yl oxc nucleons in 
the core. These nucleons are put in the lowest orbitals 
allowed by the Pauli principle; (2) from the contribution 
of the nucleons A cxc outside of the closed shell. The wave 
function describing them is obtained by diagonalizing a 



suitable effective interaction in the corresponding model 
space. It turns out that this last contribution can be 
very well approximated by allowing the A oxc nucleons to 
be placed in the lowest energy harmonic oscillator shells 
allowed by the Pauli principle. This approximation is 
quite good even in the case of light nuclei [2l| . 
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